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' Masses of hadrons (baryons as well as mesons) are modified in the nuclear medium because 

■ of their interactions. In this paper we investigate the effect of in-medium meson masses on 

the properties of the nuclear matter. The calculations are performed in the Walecka model. 
We find that with the inclusion of meson mass modification, the computed equation of state 
becomes softer. We also find that the decrease in the nucleon as well as meson masses in 
the medium is smaller than that obtained in the case when the in-medium meson masses 
are not taken into account. 
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I. INTRODUCTION 

Walecka model |j],U and its extensions, which include nonlinear meson interactions || and derivative 
coupling Qj, have been used extensively in the investigation of nuclear systems for quite some time now. 
One of the attractive features of these models is that these are relativistic models based on meson-nucleon 
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interactions. The saturation properties of the nuclear systems is in-built in these models through scalar (a) 
and vector (uj and p) mesons interacting with the nucleons. The scalar meson generates long range attraction 
and the vector mesons give rise to short range repulsion. Furthermore, the spin-orbit interaction which is 
essential in nuclear physics arises in a very natural way in these models. A number of nuclear structure 
calculations based on the extensions of the Walecka model || have been performed. These calculations have 
succeeded in explaining the properties of stable nuclei across the periodic table as well as those of nuclei far 
away from the stability line. In a nut-shell, the calculations in the Walecka model are performed in a mean 
field approximation where one solves the nucleon and meson (classical) field equations self-consistently. Thus, 
the calculations are done in Hartree approximation and Fock (exchange) terms are not included. Essentially, 
this amounts to saying that the meson fields develop nonzero vacuum expectation values in the presence of 
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nucleons and these mean fields, in turn, generate scalar and vector potentials for the nucleons. The effective 
(in- medium) nucleon mass (m* ) is then the sum of free nucleon mass and the scalar potential. 

The success of the Walecka-type models in explaining the nuclear properties has led to the application 
of these models to other areas. This essentially amounts to extrapolating the applicability of the model 
to large isospin, nuclear densities much larger than the normal nuclear matter density (pnm) and nonzero 
temperatures. For example, the model has been used to compute the equation of state of neutron matter 
which is used in neutron star calculations ||. It has also been employed in the calculations of nuclear 
equation of state at densities higher than pnm 0- Recently calculations of in- medium meson masses at 
high densities and temperatures in Walecka model have also been reported These calculations were 

motivated by the conjecture that the properties of hadrons may be modified in nuclear medium fiofl and 
this may have significant consequences on the reactions occuring in the nuclear medium. These calculations 
show that p and ui meson masses are significantly altered in the nuclear medium. For example, uj mass may 
be as low as 630MeV even at pnm which is about 0.8 times the free uj mass m w . 

The Walecka model was initially introduced as an effective model for the calculation of properties of nuclei. 
Thus the scalar and vector mesons were effective fields introduced to account for the main features of the 
N-N interaction. It is true that the pion, which generates the long range NN interaction is not included. But 
it is argued that in the mean field approximation for spin and/or isopin zero systems one pion exchange does 
not contribute and the two pion exchange is incorporated through a meson exchange. This argument breaks 
down at high temperatures where pionic excitations should dominate. If one takes this argument seriously, 
in the spirit of effective interactions of nuclear structure or the Fermi theory of weak interaction, then it is 
not meaningful to compute the in-medium meson masses in Walecka model. That is, since the parameters 
of the model (meson masses and the coupling constants) are fitted to satisfy the nuclear properties (binding 
energies, nuclear incompressibility, symmetry energy etc.), it does not make sense to compute the in medium 
meson masses again. In other words, the fitted meson masses are in medium masses and therefore should 
not be computed again. 

The argument given in the preceding paragraph is reasonable if the variation of the effective masses as 
a function of nuclear density or temperature is weak. However, the calculations || show that the effective 
masses change quite a bit as a function of density. This means that the fitted value of the meson masses at 
Pnm will change significantly at higher densities and the computation of the nuclear matter properties at 
these densities will not be reliable. 

Clearly, a better procedure would be to use the in-medium meson masses in Walecka model calculations 
particularly if the results are sensitive to the meson masses. We do believe that such self-consistent calcu- 
lations would yield results different from the standard Walecka model calculations. To illustrate this point, 
let us consider nucleons interacting with scalar mesons of mass m. The nucleon- nucleon interaction poten- 
tial generated by the scalar meson exchange is then V(r) = g 2 e ijrr where g is the meson- nucleon coupling 
constant (see for details fll|| ). If we consider the meson mass in the medium is modified to m e ff, the 
nucleon- nucleon potential in the medium is obtained by replacing m by m e f f . Clearly, this change will affect 
the computation of energy/nucleon, effective nucleon mass etc. In particular, the decrease in the meson 
mass will lead to increase in the range of the potential and decrease in its strength. Further more, since 
the meson mass in the medium depends on the nuclear density, the strength and the range of the potential 
will be density-dependent. Thus, one expects that nuclear matter properties computed using this procedure 
are likely to differ significantly from the standard Walecka model calculations particularly at densities larger 
than p NM - 
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The preceeding discussion clearly brings out the problems associated with the computation of the meson 
masses in the medium blindly. The correct procedure should be to compute meson masses as well as nuclear 
matter properties self-consistently (upto certain order in perturbation) and fit the nuclear matter properties 
by adjusting the parameters of the Walecka model. Note that this self-consistency is over and above the 
self-consistency in the standard Walecka model calculations. We have done this in the present work. Here 
we have employed the Walecka model with a and uj mesons coupling to nucleons. From our calculations we 
find that it is indeed possible to carry out such a program and fit the nuclear matter properties (binding 
energy /nucleon at nuclear matter density). We also find that the change in the nucleon and meson masses in 
the nuclear medium are less pronounced and the computed equation of state is softer. The computed nuclear 
matter incompressibility is about 460 MeV in comparison with 610 MeV, the value computed in standard 
the Walecka model. 

Before we go on to the details of our calculation we would like to note that our calculations can be repeated 
for the extensions of the Walecka model . We have not done this here since we want to focus on the 
effect of the variation of meson masses on the nuclear matter properties as a function of nuclear density. We 
believe one would obtain qualitatively similar results for the extensions of Walecka model. 

II. MESON MASSES AND NUCLEAR MATTER PROPERTIES 

The Walecka model Lagrangian, having nucleons interacting with scalar er-meson and vector w-meson is 
given by § 

C w = ^{h^ - m n )ip - g^^uj* 1 ~ ^G^G^ + ^m^ (I w' 1 

+ \(d^a - mla 2 ) + (1) 

In the above equation tp, a and w M are, respectively, the nucleon, the scalar and the vector meson fields; 
m^m^ and m w are the corresponding masses; g a and g^ are the couplings of the nucleon to scalar and 
vector mesons, respectively and G M „ = d^uj v — d v uj^. 
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FIG. 1. Diagrams contributing to nucleon self energy ( a ) and a and ui self energies ( b ). 

In the following discussion we have described the study of the in-medium properties of hadrons using the 
above Lagrangian. The in medium nucleon mass, the energy density per baryon and the a and uj masses are 
calculated self-consistently. The coupling constants are determined by fiting the nuclear matter saturation 
density to be pnm = 0.15/m~ 3 , the corresponding binding energy being —15.75MeV per baryon. 
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The effective nucleon mass in the Relativistic Hatree Approximation (RHA) can be calculated from figure 
la. This is given by || 
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where G is the in-medium baryon propagator which is given by 

1 



G(p) = (7^+<) 



v 



2iirS(p 2 - m* n z )6( Po - k F ) 



(2) 



(3) 



So we get 



rn, 



9l 



*3 „,3\ 



k F Ep + rri^ln 
7 



k F + E F Y 


9l 


mfZn ( 




m n J. 


7r 2 m* 2 




13 

- —m n (m 


* \2 i 


25 

12 (m " 




-"In; 



(4) 



The effective mass of a meson can be calculated from the solution of the Dyson-Schwinger equation which 
relates the free and the full propagator as follows : 



D(p) = D (p) + D (p)UD(p) 



(5) 



where Do(p) is the free Green's function, D(p) is the full Green's function and II is the polarisation function 
(at zero three momentum). The pole of the full propagator then gives the effective mass of the meson. 

We calculate the meson masses in Random Phase Approximation (RPA). The corresponding diagram is 
given in figure lb. The effective mass of the a meson is given by [121 
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The polarisation function for the cr-meson has two parts. First part is the in-medium contribution and 
the second part is the vacuum contribution. In calculating the vacuum part we have renormalised the 
contribution according to ref. Jl3| . 

For the omega meson we have 
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where II^t is the transverse part of the polarisation tensor (in the static limit the transverse and the 
longitudinal components are the same) which is given by 
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The transverse part is given by 
IW = 
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The binding energy per baryon, calculated in RHA, is given by, 
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where £ is the energy density which is (T 00 ) and T^ v is the energy momentum tensor. The energy density 
in the RHA is given by 
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We now solve the above equations self-consistently. The parameter values that we get, from the nuclear 
matter saturation properties, are given by g a = 6.716 and g u = 7.64. In case we do not consider the 
in-medium meson masses the corresponding values are 8.69 and 9.9 respectively. The results are shown in 
figures 2-3. In figure 2 we have plotted the masses of nucleon, scalar meson and vector meson as a function 
of density. In order to make a comparison, both the cases i.e. with and without self-consistency have been 
plotted. In figure 3 the energy per baryon has been plotted for both cases. 



III. RESULTS AND DISCUSSION 



We now come to the discussion of the results. Figure 2 displays the plots of hadron masses as a function 
of nuclear density. The dashed lines correspond to the standard Walecka model calculation where the 
modification of meson masses are not included in the calculation of nuclear matter properties and the 
continuous lines correspond to our (self-consistent) calculation. We find that in a self-consistent calculation 
the change in the masses is somewhat less dramatic. In particular, the change in the in-medium a and 
uj masses is less than 15% in a self-consistent calculation whereas the maximum change is more than 20% 
for the standard Walecka model. Hence at nuclear matter density the a and u> masses are 471MeV and 
686MeV respectively for the self-consistent case, whereas the corresponding values for the case where there is 
no self-consistency are 445MeV and 632MeV respectively. This in terms of the N — N interaction potential 
implies that the potential is more attarctive for the self-consistent case. As a result the nuclear matter will 
be more compressible i.e. the incompressibility will be lower. Such changes, in general, will have bearing 
on the binding energy as well as on the equation of state of the system. The behavior of the nucleon mass 
appears to be similar for both the cases but the effective nucleon mass is slightly larger in a self-consistent 
calculation ( 0.79 tojv vs 0.73 tojv at pnm )■ From the behavior of the masses one can conclude that in order 
to obtain a reliable estimate of masses of hadrons in nuclear medium it is essential to do a self-consistent 
calculation. 
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FIG. 2. Nucleon, ui and a masses effective masses. Continuous and dashed lines are for with and without 
self-consistancy. 

Figure 3 shows a plot of the binding energy/nucleon ( B/A ) as a function of the nuclear density. Both the 
curves have a minimum at the nuclear matter density, simply because the binding energy at the saturation 
density is fitted (— 15.75MeV at 0.15.fm~ 3 ). One should, however, note that the second derivative of B/A, 
which is related to the nuclear incompressibility, is smaller for the self-consistent calculation. The computed 
incompressibility for the self-consistent calculation is 460 MeV in comparison with the standard Walecka 
model result of 610 MeV. In other words, the self-consistent calculation drives the incompressibility towards 
lower values. The binding energy per baryon is much lower at higher densities for the self-consistent case. 
This implies that the equation of state obtained in the self-consistent calculation is softer. This is an 
important result in the context of neutron star physics. In the standard Walecka model the equation of state 
is stiffer resulting a value of the neutron star mass which is higher than the experimental bound. Our results 
indicate that a self-consistent calculation may yield a better estimate for the neutron star masses. 

So, to conclude, we have calculated the nuclear matter properties as well as the in-medium meson masses 
self-consistently by including the change in the meson masses due to medium effect in the calculation of 
nuclear matter properties. Our study gives rise to two important conclusions. First of all, we find that 
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the change in the meson masses due to the effect of the medium is smaller in a self-consistent calculation. 
This means that a calculation which does not impose the self-consistency over estimates the meson masses. 
The second result is that the self-consistent calculation yields a softer equation of state and lower values of 
incompressibility. We believe that this qualitative behaviour will persist for asymmetric nuclear matter and 
will have a bearing on neutron star properties. 
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FIG. 3. Binding energy /nucleon ( B/A ) as a function of nuclear density. The continuous and dashed curves are 
for with and without self-consistancy 

We would like to emphasize again that, even if the meson fields are considered as effective fields, the 
self-consistency described in the present work is necessary if one wants to apply the Walecka model at higher 
densities. On the other hand if these fields are considered as real fields, it is essential to do a self-consistent 
calculation to obtain a reliable estimate of meson masses in the medium. It would be interesting to investigate 
the behaviour of p meson mass in the medium since it has a bearing on the quark-gluon plasma diagnostics. 
Such a calculation is in progress. 
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